Quantum particle escape from a time-dependent confining potential 
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Quantum escape of a particle via a time-dependent confining potential in a semi-infinite one-dimensional 
space is discussed. We describe the time-evolution of escape states in terms of scattering states of the quantum 
open system, and calculate the probability P(t) for a particle to remain in the confined region at time t in the 
case of a delta-function potential with a time-oscillating magnitude. The probability P(t) decays exponentially 
in time at early times, then decays as a power later, along with a time-oscillation in itself. We show that a 
larger time-oscillation amplitude of the confining potential leads to a faster exponential decay of the probability 
P(t), while it can rather enhance the probability P(t) decaying as a power. These contrastive behaviors of the 
probability P(t) in different types of decay are discussed quantitatively by using the decay time and the power 
decay magnitude of the probability P(t). 

PACS numbers: 05.60.Gg, 03.65.Nk, 73.23.-b 



I. INTRODUCTION 

Quantum escape, as a quantum mechanical leakage of ma- 
terials from a confined area, is an important concept to study 
dynamical properties of quantum open systems. It appears in 
many physical phenomena, such as the a decay of a nucleus 
IHI2I, radiative decay of molecules fU], transition of states in 
chemical reactions 0], etc. Studies on quantum escapes have 
been done in one-dimensional systems with a stationary po- 
tential barrier &0L billiard systems with leads H], kicked 
rotator models Istil lll. n etwork systems Il2l4l4tl . and many- 
particle systems Il5l - [l7ll . and so on. 

In escape phenomena, materials continue escaping from a 
confined area as an irreversible process, so various quantities 
decay in time. As one of such decaying quantities, the prob- 
ability for particles to remain in the confined area, which we 
call the survival probability d [3 H [H [H in this paper, 
has been widely investigated in escape phenomena. In quan- 
tum escapes, the survival probability decays exponentially 
in time 

U Gj, H HI, by tunneling via a potential barrier, 
etc. Moreover, it also shows a power decay after a long time 
dd[3,|2lj,|22l. As decay properties of the survival probabil- 
ity, time scales of its exponential decay l9t 4Tlll . and changes of 
the power of decay under different conditions |[T^[l6l[T7ll23ll . 
etc. have been studied. 

The principal aim of this paper is to investigate effects of 
a time-dependent confining potential in quantum escape phe- 
nomena. Quantum dynamics of systems with a time-periodic 
potential is analyzed by the Floquet theory, etc. Ir24ll25tl . and 
recent experimental and theoretical developments in micro- 
scopic or mesoscopic systems have pushed to study quantum 
dynamical behaviors caused by a time-dependent external ma- 
nipulation l2oT-|}28tl . In this paper, we consider quantum es- 
capes of a particle via a time-periodic confining potential in 
a semi-infinite one-dimensional space, as a model of a time- 
dependent manipulation of particle escapes. In the case of 
a time-dependent delta-function potential, we obtain an an- 
alytical form of the quantum propagator to describe a time- 
evolution of the wave function, leading to a calculation of the 
survival probability of the system at any time. We show that 
a time-oscillation of the confining potential has opposite ef- 



fects in different decay behaviors of the survival probability. 
In the exponential decay of the survival probability, a larger 
time-oscillation amplitude of the confining potential leads to a 
faster decay of the survival probability. In contrast, the time- 
oscillation of the confining potential can rather enhance the 
survival probability decaying as a power. We calculate the de- 
cay rate of the survival probability and the magnitude of its 
power decay, by which these reduction and enhancement of 
the survival probability by a time-oscillation of the confining 
potential are discussed quantitatively. 



II. QUANTUM OPEN SYSTEMS WITH A 
TIME-DEPENDENT POTENTIAL 

A. Semi-infinite one-dimensional systems with a 
time-periodical potential 

The system, which we consider in this paper, is a quantum 
particle system with a time-dependent potential in the semi- 
infinite one-dimensional space. The wave function ^(x, t) of 
the system at position x and time t satisfies the Schriidingier 
equation ihd^(x,t)/dt = H(t)&(x,t). Here, H(t) = 
— [ti 2 /(2m)]d 2 /dx 2 + U(x,t) is the Hamiltonian operator 
with the mass m of particle, Planck's constant 2irh, and the 
potential U(x,t) at position x and time t. We impose the 
hard-wall boundary condition at the end x = of the semi- 
infinite one-dimensional space x > 0, so that the wave func- 
tion &(x, t) satisfies the condition ^(0, t) = at any time t. 
We assume that the potential U(x, t) is a periodic function of 
time, so that 



+00 



U(x,t) = U(x,t + T)= 



(x)e i 



(1) 



n— — 00 



with the time period T = 2-k/lu, in which u is a posi- 
tive constant and u^ n \x) is a function of position x only 
in = 0, ±1, ±2, • • ■)■ We impose the condition vS~ n '{x) — 
[uW(i)| for the function of x in Eq. |[T), so that 

the potential U(x, t) is a real function of x and t. Here, X* 



2 



with the asterisk * denotes the complex conjugate of X for 
any complex number X. 

We consider the scattering state $k(x, t) at position x and 
time t, which is induced by the incident plain wave from the 
infinite region x — > +00 with the energy Ek = (hk) 2 /(2m) 
and the wave number k. As a solution of the Schriidingier 
equation with the time-dependent potential ([T), the scattering 
state is represented as 



with T^™' defined by 



+ OC 



<t> { k\x)e~ l(Ek+nh " )t,h , 



(2) 



T< n) = k 2 



2mu> 



(4) 



where the function </> l"^ (x) o f x satisfies the time-independent 
differential equation II241I26I1 as 



d 2 ^\x) 



dx 2 



2m 
J 2 



E < 

f-i— OO 



(3) From Eq. © we derive the Lippmann-Schwinger equation 



+00 



4*° 0) = S n0 \j z sin{kx) + ^ Y 



dy G^\x,y)u^- n \y)^\y) 



(5) 



where { x i y) is given by 



B. Scattering states with a time-dependent delta-function 
potential 



G*?\x,y) 



>\x-y\ 



'\x+y\ 



2h r 



as a Green function satisfying the differential equation 



d 2 G[ n \x,y) 
dx 2 



T^GP(x,y)=5(x-y) 



(6) 



(7) 



From now on we especially consider the case of 

u W(x) = \ in) 8(x~L) (8) 

with constants L (> 0) and \( n \ satisfying the condition 
A(~") = [AW]*, n = 0,±1,±2,---. Under the condi- 
tion ((HJ the potential is represented as U (x, t) = d(x — 
L ) Y,n=-oo X^e inut , i.e. the delta-function potential with 
a time-dependent magnitude. By Eq. (O, Eq. @ is trans- 
formed into 



with the boundary condition g£ (0,y) = 



0. Here, in the scattering state (f2]i we imposed the con- 
dition that if there is no potential energy then the scatter- 
ing state 5>fc(a;, t) is represented as the state ip^Xjt) = 
^2/irsm(kx) exp(—iEkt/h) induced by the incident plain 
wave with the wave number k, with the orthogonal relation 
fo°°dx [i()k{x,t)]*il)k>( x ,t) = 5(k - k') for kk' > and 
the boundary condo ipk(0,t)=Q. This condition determines 
the form <5„o-\/2/7r sin(fca;) of the first term in the right-hand 
side of Eq. (|5}. We also defined the square root of Y^, as 

^/Tp = iv/|Tjp| for the case of t[' i) < in Eq. ©. 
More explanations of Eq. (O and the Green function (|6]l are 
given in Appendix lAl 



'Pk"^ ( x ) = $no\ — sin(fcx) 

V 7T 



~G^(x,L) £ ^- n) ^(L) (9) 



fi— —00 



without any spatial integration. 

We in 
element 



We introduce the matrix Aj, = (Al ) with the jiv matrix 



Ai"^^G^(L,L)A(-"). 



(10) 



Using the matrix element ([Tol l, Eq. (O at x = L leads to the 



relation Y^t=- 



^\L) = 5 n0 J^sm(kL), 
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namely 



(L) = J-sm(kL)[(I-A k ) 



-11 (nO) 



(ID 



with the identical matrix /. Here and hereafter we put X^"' 
as the fii/ element of X for any matrix X in this paper. By 
inserting Eq. ( TTTT i into Eq. (O and using Eq. ( fTOb we obtain 
the quantity <p£> (x) as 



4>k^(. x ) = 6no\l - sin(kx) 



* G ( k\L,L) 



-li(nO) 



(12) 



The scattering state 5>fc(x, t) induced by the incident plain 
wave with the wave number k is given by Eq. (O with Eq. 

It is essential to note for contents of this paper that the scat- 
tering state $fe(x, t) satisfies the orthogonal relation as 

r + CO 

/ dx [$ k (x,t)]*$ k ,{x,t)=6(k-k') (13) 
Jo 

for any time t and almost all positive numbers k and k . The 
proof and explanations of Eq. ( fT3l are given in Appendix IB1 



C. Time-evolution of the wave function with the propagator 

In this paper we restrict ourselves our considerations into 
the quantum state *ff(x, i) which is expanded by the scattering 
state 3>fc(x, t) as 

/•+oo 

*(*,*) = / dkA k $ k {x,t) (14) 
Jo 

with the coefficient A k . By using Eqs. ( fT3l and ( TBI , the inner 
product J Q + °°dx [$ fc (x, 0)]* *(x, 0) of 3> fe (x, 0) and *(x, 0) 
leads to the coefficient namely 



dx [$ fc (x,0)]**(x,0). 



(15) 



By inserting Eq, ( fT3T > into Eq. ( fT~4T >. the wave function 4'(x, f) 
at position x and time t is rewritten as 

r + oo 

$(x,t) = / dy K(x,y;t)V(y,0), (16) 
Jo 

where if(x, y; t) is defined by 

K(x,y;t)= / dfc$ fc (i/,0)*$ fc (i,t) (17) 
Jo 

as the propagator for the time-evolution of the wave function 
*(x,t). 



III. PARTICLE ESCAPE FROM A TIME-OSCILLATING 
CONFINING POTENTIAL 

A. Survival probability 

Now we consider the initial condition in which the particle 
exists only in the finite region (0, L) at the initial time t = 0, 
namely 



*(x,0) 



* (x) for < x < L 
for L < x 



(18) 



with a function ^q(x) of x. Staring from this initial condition, 
we investigate the probability P(t) for a particle in the finite 
region (0, L) at time t as 



P{t) = [ dx \^{x,t)f 
Jo 



(19) 



which we call the survival probability. 

Furthermore, from now on we assume that the initial state 
is given by the energy eigenstate of the particle confined in the 
finite region [0, L], namely, 



*o(a0 



2 /07T 

L sm \T x 



(20) 



for < x < L with an integer value a = 1, 2, • • • . Under 
this specification of the initial conditions ( fT8l and d20b we can 
carry out analytically the integral over the position x in Eq. 
( TT9b . We give an explicit form of the survival probability P(t) 
in Appendixlclunder this initial condition. 

As a concrete time-dependence in the potential U(x,t), in 
the next subsection llll Bl we consider a time-oscillating poten- 
tial, i.e. 

U(x,t) = [a + bcos(u>t)]5(x- L) (21) 
with real constants a and b. This potential is derived from 



a for n = 
| for n = ±1 
otherwise 



(22) 



for the constants A*-") to specify the potential U(x, t) by Eqs. 
© and ©. 



B. Time-oscillating potential effects in decay of the survival 
probability 

In this subsection we show graphs of the survival probabil- 
ities P(t) as functions of time t under the initial conditions 
( fT8l and (|20l > with the time-dependent potential ( f2Tb for vari- 
ous values of the parameters a and b. For calculations of such 
survival probabilities we use the parameter values as L = 1, 
uj = 2tt [corresponding to the time -period T = 2ir/u) = 1 
of the potential (fJTTil. a = 1, m = 1 and H = 1. In actual 
numerical calculations we replaced the integral with respect 
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FIG. 1: (Color online) The survival probabilities P(t) as func- 
tions of time t for the cases of (a, b) = (1.5, 2) (the solid red line), 
(a, b) = (1.5, 1) (the dashed green line) and (a, b) = (1.5, 0) (the 
dotted blue line). The main figure is the linear-log plots of P(t), 
and the inset is their log-log plots as well as the asymptotic power 
decay £/t 3 corresponding to the case of (a, b) = (1.5, 0) shown in 
the straight shin black line. Here, and in all figures hereafter, we use 
dimensionless units with L = 1, m = 1 and h = 1, 



FIG. 2: (Color online) The survival probabilities P(t) as functions 
of time t for the cases of (a, b) = (5, 4) (the solid red line), (a, b) = 
(5, 2) (the dashed green line), (a, b) — (5, 0) (the dotted blue line), 
(a, b) = (5, —2) (the dash-dotted purple line) and (a, b) — (5, —4) 
(the dash-double-dotted cyan line). The main figure is the linear-log 
plots of P(t), and the inset is their linear-linear plots up to the first 
oscillating period T = 2n /to = 1 of the potential. 



to the wave number k up to +oo by the integral with respect 
to the wave number k up to 20, checking that the absolute 
value of the actual integrand is extremely small at k = 20 
so that the contribution of the integral over k > 20 would be 
negligible. For numerical calculations of the probability P{t) 
we also have to replace the infinity times infinity matrix A& 
with a finite (27V + 1) times (27V + 1) matrix (A^) with 
// = -TV, -7V+ 1, ••■ ,7V and v = -TV, -7V+ 1, • ■ ■ ,7V. 
This implies that we take into account of effects up to the 7V- 
th closest energy scattering states from the state with the wave 
number k for a small magnitude of the time-depending part 
of the potential ( f2TT >. Numerical calculations of the proba- 
bility P(t) with a large dimensional matrix (A! ) is very 

time-consuming, and we used the 7 times 7 matrix (A^-*), 
i.e. TV = 3, for the numerical results in this subsection 1111 B I 
For some cases discussed in this subsection we calculated the 
probabilities P(t) not only for TV = 3 and but also for TV = 4 
with the same values of the other parameters. Differences be- 
tween these calculated results for TV = 3 and TV = 4 were al- 
most negligible comparing to values of the probabilities P(t) 
themselves, although these results were partly different from 
the ones for TV = 2. These facts would justify to calculate the 
probability P(t) by TV = 3 for values of the parameters a and 
b, etc., shown in this subsection HIIBI 

Figure Q] is the survival probabilities P(t) as functions of 
time t for the cases of (a, b) = (1.5, 2) (the solid red line), 
(a, b) = (1.5, 1) (the dashed green line) and (a, b) = (1.5, 0) 
(the dotted blue line), as their linear -log plots (the main figure) 
and their log-log plots (the inset) 12911 . Here, in its inset we 
added the shin black line showing the asymptotic power decay 



of the survival probability in time for the time-independent 
potential with b = 0, as given analytically by P(f) L_ n ^'i 00 
i/t 3 with i ee 4,m 3 L 6 h 5 /[3<T 2 Tr 3 {h 2 + 2mLa) 4 } UM. It is 
shown in this figure that the survival probability P(t) itself 
oscillates in time for the case of b ^ by a time-oscillation 
of the potential ( f2TT >. Figure Q] also shows that there are two 
types of decay behaviors of the survival probability P(t); one 
is an exponential decay at early times, and the other is a power 
decay later, apart from its time-oscillatory behavior. More- 
over, a comparison of the time-independent potential case 
of (a,b) = (1.5,0) with the cases of (a,b) = (1.5,1) and 
(1.5, 2) in Fig. Q]suggests that the time-oscillation of the po- 
tential reduces values of the survival probability P(t) decay- 
ing exponentially, while it can increase values of the proba- 
bility P(t) decaying as a power. In the following subsections 
III1 B II and HUB 21 we will discuss these behaviors caused by 
a time-oscillating potential in different types of decay of the 
survival probability separately in more details. 



1. Exponential decay of the survival probability 

Figure [2] is the survival probabilities P(t) as functions of 
time t for the cases of (a, b) = (5, 4) (the solid red line), 
(a, 6) = (5,2) (the dashed green line), (a,b) = (5,0) (the 
dotted blue line), (a, b) = (5, —2) (the dash-dotted purple 
line) and (a, b) = (5, —4) (the dash-double-dotted cyan line). 
Here, the main figure is the linear-log plots of P(t), and the 
inset is their linear-linear plots at beginning times up to the 
first oscillating period T = 2-k/lu = 1 of the potential d2ll . 
For this figure, we chose a larger value of the parameter a 
than that in Fig. Q] so that the survival probabilities show their 
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FIG. 3: The decay time r of the survival probability decaying expo- 
nentially as a function of the parameter b of the potential U(x, t) = 
[a + bcos(uit)]S(x — L) for the case of a = 5. The dashed line is 
a fitting quadratic curve r = C,i — (,2b 2 for the numerical data of r 
with fitting parameters £1 and £2- 



FIG. 4: (Color online) The survival probabilities P(t) as func- 
tions of time t for the cases of (a, b) = (0.5, 2) (the solid red line), 
(a, b) = (0.5, 1.2) (the dashed green line), and (a, b) = (0.5, 0) (the 
dotted blue line), (a, b) — (0.5, —1.2) (the dash-dotted purple line) 
and (a, 6) = (0.5, —2) (the dash-double-dotted cyan line) on log-log 
plots. The straight shin black line shows the asymptotic power decay 
£/t 3 for the case of b = 0. 



exponential decays more clearly. 

The almost straight lines in the main figure of Fig. [2] as 
linear-log plots of P(t) show that at the presented times the 
survival probabilities P(t) decay exponentially in time, apart 
from their small time oscillations with the period T. These ex- 
ponential decays continue for a longer time for a larger value 
of the parameter a of the potential (|2TT >. suggesting that it is 
caused by a quantum tunneling via the delta-function poten- 
tial. The survival probability P(t) for 6^0 oscillates in time 
with a larger amplitude for a larger time-oscillation amplitude 
|6| of the potential (f2TT >. The survival probabilities with the 
same absolute values of the parameter b, such as b = ±2 or 
b = ±4, decay exponentially in time with the same decay rate 
but with the opposite phase of their time-oscillations. 

The time-dependent potential barrier magnitude a + 
bcos(ujt) of the potential d2TT > is larger for a larger value of 
b for < t < T/4 (and for 3T/4 < t < T), so that the sur- 
vival probability P(t) decays slower for a larger value of b, as 
shown in the inset of Fig. |2] After it, for T/4 < t < 3T/4 the 
potential barrier magnitude becomes to be smaller for a larger 
value of b, so that the survival probability P(t) should decay 
faster for a larger value of b at this time. As shown in the inset 
of Fig. |2 the magnitude of such an enhancement of the decay 
speed of the survival probability P(t) in T/4 < t < 3T/4 
is larger than that of its suppression in < t < T/4 and 
3T/4 < t < T for a larger value of b. This kind of unbalance 
of decay speeds of the survival probability P(t) for each os- 
cillating period T of the potential (|2TT > causes a faster decay of 
the probability P(t) on average for a larger value of the time- 
oscillation amplitude \b\ of the potential d2"TT i as shown in the 
main figure of Fig. [2] 

In order to discuss time-oscillating effects of the potential 
(I21T 1 in exponential decays of the survival probabilities, we fit 
each exponential decay of the survival probability P(t) to an 



exponential function 7cxp(— t/r) with the fitting parameters 
7 and t, and introduce the decay time r as its fitting value. 
Figure[3]is a graph of such a decay time r as a function of the 
parameter b for the case of a = 5. This figure suggests that the 
decay time t is a decreasing function of the time-oscillation 
amplitude \ b\ of the potential, and is also an even function of b. 
As shown in the dashed line of this figure, the 6-dependence 
of the decay time r is well fitted to a quadratic curve r = 
Ci - C2& 2 with constants Ci = 2.61 and C2 = 8.11 x 1CT 2 
near 6 = 0, although a deviation from it is recognized for a 
large value of |6|. 



2. Power decay of the survival probability 

In Fig. |4]we show the graphs of the survival probabilities 
P(t) as functions of time t for the cases of (a, b) = (0.5, 2) 
(the solid red line), (a, b) = (0.5, 1.2) (the dashed green line), 
(0,6) = (0.5,0) (the dotted blue line), (a, b) = (0.5,-1.2) 
(the dash-dotted purple line) and (a, 6) = (0.5,-2) (the dash- 
double-dotted cyan line). This figure is their log-log plots 
to show their power decays as straight lines, and also shows 
the asymptotic power decay £/i 3 obtained analytically for the 
case of b = 0. Here, we chose a smaller value of the parame- 
ter a than that in Fig. [1] so that power decay behaviors of the 
survival probabilities appear more clearly. 

As shown in Fig. |4j after a finite time (approximately t > 5) 
the survival probabilities P(t) decay as a power oc i~ 3 , apart 
from their time-oscillations. The survival probability in such 
a power decay has a tendency to take a large value for a large 
value of the time-oscillation amplitude |6| of the potential 
d2Tb . but such increasing rate of the survival probability for 
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FIG. 5: The power decay magnitude r\ of the survival probability 
decaying as a power as a function of the parameter b of the poten- 
tial U(x,t) = [a + bcos(u)t)]5(x — L) for the case of a = 0.5 
on a linear-log plot. The dashed line is a fitting curve given by 
11 = ?i exp [ft (6 — ft) 2 ] for the numerical data of 77 with fitting 
parameters ft , ft and ft . 



< is different from that for b > 0. The time-oscillation 
of the potential also causes a time-oscillation of the survival 
probability P(t), whose amplitude is larger for a larger value 
of I fo| . The phase of time-oscillations of the survival probabil- 
ities for b < is opposite to that for b > 0. It may be noted 
that the time-oscillations of the survival probabilities in their 
power decays involve clear temporal increases of the value of 
P(t) as the time t passes, differently from the ones in their 
exponential decays. 

In order to investigate quantitatively enhancement of the 
survival probabilities by the time-oscillation of the potential 
(l2lT i in their power decays, we fit them to the power decay 
function ?/<~ 3 with the fitting parameter 77 as the power decay 
magnitude, apart from their time-oscillations. Figure [5] is a 
graph of such a fitting value of 77 as a function of the parame- 
ter b of the potential (1211) . Here, we used the analytical value 
77 = £ w 2.69 x 10~ 3 only for the case of b = 0. This fig- 
ure clearly demonstrates a tendency in which the power decay 
magnitude 77 of the survival probability decaying as a power 
takes a large value for a large time-oscillation amplitude \b\ 
of the potential ( f2TT >. However, it is important to note that the 
power decay magnitudes 77 is not an even function of b, and it 
takes the minimum value at a non-zero negative value of b. To 
clarify this point, we fitted the numerical data of 77 around its 
minimum point to the function r) = ft exp [ft (6 — <s) 2 ] with 
fitting parameters ft, ft and ft. Here, we used a quadratic fit- 
ting function to the data of log 77 rather than of 77 itself for a 
better fit for a wider values of b, and the used values of the fit- 
ting parameters are ft = 2.20 x 10~ 3 , ft = 6.36 x 10 _1 and 
ft = —5.09 x 10 _1 . This function fits our numerical results 
of 77 reasonably well around its minimum point, and suggests 
that the power decay magnitudes 77 takes its minimum value 
at b « —0.5 approximately, for the case of a = 0.5. 



IV. CONCLUSION AND REMARKS 

In this paper we discussed quantum escapes of a particle 
via a confining potential with a time-dependent magnitude in 
a semi-infinite one-dimensional space. We calculated scatter- 
ing states of the quantum open system, and based on their or- 
thogonal relations we derived the quantum propagator for the 
time-evolution of wave function of the system with a time- 
dependent potential. Our special interest in such a quantum 
open system was decay properties of the survival probability, 
i.e. the probability for a particle to remain in the finite re- 
gion between the wall and the confining potential. By using 
the propagator represented by the scattering states, we calcu- 
lated the survival probability for the case of a delta-function 
confining potential with a time-oscillating magnitude, and the 
initial state as the energy eigenstate of a particle confined in 
the finite region. The survival probability shows two different 
types of decay behaviors: an exponential decay at early times 
and a power decay later. The time-oscillation of the potential 
magnitude reduces a decay time of the survival probability de- 
caying exponentially. In contrast, a large time-oscillation of 
the confining potential enhances a magnitude of the survival 
probability decaying as a power. The power decay magnitude 
of the survival probability is not symmetric under change of 
the sign of the time-oscillating part of the potential, different 
from the decay time of the survival probability decaying ex- 
ponentially. 

In this paper, we calculated numerically the double inte- 
grals over wave numbers for the survival probability. It may 
be noted that such numerical integrations for a time-dependent 
potential are much harder than the one for the corresponding 
time-independent potential because of more complicate time- 
dependence of the propagator and matrix treatments. This 
kind of difficulties required us some additional numerical 
techniques (such as a calculation of the matrix Ap(I — A^) -1 
as solutions of linear algebraic equations, etc.), and actually 
restricted us parameter values under which we could calcu- 
late the survival probability within a reasonable calculation 
time and a numerical accuracy. It would be a future problem 
to improve further numerical calculation methods, as well as 
analytical ones, to calculate the survival probability for wide 
range of values of various parameters. 
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Appendix A: Lippmann-Schwinger equation in the semi-infinite 
one-dimensional space 

In this appendix, we explain the Lippmann-Schwinger 
equation (O with the outgoing Green function (|6]i as the solu- 
tion of the differential equation ©. 



7 



We note that the function sin(fca;) of k and x(> 0), pro- 
portional to the energy eigenstate with a wave number k at 
position x for the case of no potential in the semi-infinite one- 
dimensional space (0, +oo), satisfies the relation 

1 f + °° 

S(x — y) = — / dk sin(fcx) sin(fcy) (Al) 

* J-oo 

I 



for any real numbers x and y satisfying the condition xy > 0. 
By Eq. (lAlb the Green function (x, y) satisfying Eq. (0, 
corresponding to outgoing scattering states, is given by 



T 



(n) 



dk sin(fcx) sin(fcy) 



^dx 2 

1 r+oo ik'(x+y) , p ik'{x+y) _ ik' (x—y) _ p —ik' (x—y) 

— lim / dk' 

47T e->+0 



k' 2 



,(n) 



(A2) 



We further note the integral formula 



lim 



T (n) _ le 



,(n) 



(A3) 



where we define the square root of T 

,(n) 



(n) 



(») 



T 



TL | for the case of Tj. < 0. By carrying out the inte- 
gral in Eq. (IA21 > using the formula (IA3b we obtain Eq. ||6). 
For any function £>o(^) satisfying the equation 



d 2 g Q (x)/dx 2 + r^>g Q (x) = 0, the function <j)£'{x) = 
Qo(x) + wT,i=-°oJo°° d y Gt\x,y)u^{y)^\y) 



(«)/ 



satisfies Eq. (0). By taking qq(x) = S n o y \ sin(fcx) in this 
function ^ [x] as a special case, we obtain Eq. (0. 



for any complex number k whose imaginary part is non- 
negative. Here, d(n) and x( K ) are defined by 



d(n) 



lim — ■ 



>+o it K 2 + e 



(B2) 



x(k) = lim 



(B3) 

+0 K z + (~ ' 

is an even function of k, while 



u r " 2 J- c 2 

respectively. The function 
the function ( IB3b is an odd function of n. It is also noted that 
the function d(/t) becomes the delta function = S(k) for 
any real number n = k. We also note 



dx e 1 



AkL\ 



(B4) 



Appendix B: Orthogonal relation of scattering states with a 
time-dependent potential 

In this appendix, we give a derivation of Eq. dT3l on the 
orthogonalization of the scattering state $fe(x, t) of the system 
with the time-dependent potential given by Eqs. (Hi and (8j. 

As preparations of later calculations we note the mathemat- 
ical identity 

+oo 



(Bl) 



for any complex number k. 
We note the relation 



dx sin(ax) sin(/3a;) = — [d(a — 0) — d(a + /?)] 



ird{K) + «x(k). 



(B5) 



for any real numbers a and (3. Moreover, by using Eqs. (1BU 
and (1B41 > we can show 
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dx sin(ax)Gp (x, L) 



d\ a 



T 



- d 



- sin(aL) W [a, JT 



(B6) 



dx 



'P 



GM(L,L)V '(&\L,L)d 



^(m) 

J- a. 



Itf 



W 



T^]-Gf{L,L)w(Jr { ;\ 



(B7) 



r 



Here, we used the function W(X, Y) of X and Y defined by using the relation x(A —Y) = — x(Y — X). 

X(X-Y)- X (X + Y) 



W(X, Y) 

so that we obtain 

x(X-Y)+x(X + Y) 
2X 



2Y 



-W{Y,X) 



(B8) 



By Eq. © the inner product J Q °° dx [$> a (x, t)]* $^(x, t) 
(B9) of the scattering states is represented as 



+00 



+00 +00 



dx [$ a (x,t)]* Qpfat) = J2 J2 e i[Ea-E +(^ 



-i/)Huj]t/n I ^ x 



fji— — 00 U——00 



(BIO) 



Further, by Eqs. (fill) . ( lB5i dB6l i and (IbTT i the inner product J + °° dx cj)^ (x) <j>^ (x) in the right-hand side of Eq. (IB 101 ) is 
given by 



+00 



dx 



2 ;- + oo 

^juO^i/O - / sin(ax) sin(/3x) 



[A^J-A^)- 1 ]^ sin(/3L) 2 



+4 



+00 



dx sin(a.x)G a (x, L) 



7T JO 



[A Q (7 — Aq)^ 1 ] sin(ai) 1 2 



G^(L,L) 
1 Oo) , 



7T 



| f [A a (I-A a )-^ m S in(aL) \ [A P (I - A p )^f Q) sinQSL) 2 ^ + °° rf 

1 gL m) (l,l) J g { ;\l,l) 



G^{x,L)X G { ;\x,L) 



S^oSuO [d(a -13)- d(a + /?)] + sin(aL) sm(/3L) {A ( ( ^\ + 



(Bll) 
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where and B^"l are defined by 



A 



B 



2e 



+iS, 



fiO' 



sm(aL) 



d\ a 



d a + \ T 



-id 



i/O " 



sin(/3L) 



d /3 



-OO 



K0) ~ 



^(J-Ag)- 1 ] 

G^\L,L) 



( m o) 



[A^J-A^^WlVTfMVT 



-(m) 



[A^J-A^)- 1 ] 



_n("0) 



W a,i/T 



■/3 



[A^I-A^)- 1 ] 
G^\L,L) 



( M o) 



"00 



0, ( x/T^ 



(B12) 



(B13) 



Here, the quantity ^-{a% includes the type of for any function X(k) of complex number k because of d(a — 



T^j ±y/T^'j and the quantity in- 
cludes the terms with the function W(x, y). 

Now, we note that for the function ( lB2b we have 



?(m,"0 



0) = for any a ^ 0. Noting Eq. ( IB 14b we obtain 



X(a)d(a -/3)= X(/3)d(a - f3) 



(B14) 



G%\L,L)d(a±J?%>)=± 



G^(L,L) ' G [ ;\L,L)d[ ( \/T^ ) - Jr 



sin(aL) 



,(!/) 



d a± \/T 



-00 



G^L) 



-00 



2n/T 



,(!/) 

"(8 



2* ( vtL m) 



-00 



-a/T 



■i { [<#> (L, L)J - G^ (L, L)) D ( ( y/?M 



Here, we used the relation 



{[g£Xl,l)]*-g%\l,l)} 



-j-00 



,(!/) 



= 



(B15) 



(B16) 



(B17) 



(B18) 



to derive Eq. ( IB 17b from Eq. ( IB 16b by Eqs. © and ( IB 14b , and also the function D(X, Y) of X and Y defined by 

d(X -Y) +d(X + Y) 



D{X,Y) 



2Y 



(B19) 
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and for Eq. ( IB 15b we adopt the convention that any equation 
containing the symbols ± on the left- and right-hand sides, de- 
note two equations, one with only the upper symbol (+ in ±) 
and the other with only the lower symbol (— in ±). Moreover, 
we note the relation 



By Eqs. dBBT l, (iBTTl and (|B20T > the quantity dBHI i is rewrit- 



ten as 



[A a (I A a ) _1 ] ( ^ 0) = -S,o + [(I A*)" 1 ] [pv > .(B20) 



n (mo) 



A 



-2i 



{[(I - A a )-^} 



G%\L,L) 



(uO) 



[A^I-Ag)- 1 ] 

G ( £\L,L) 



( m o) 



[{I-A p ) 



+2^(5^0 



[Aq,(7 — Aq,) -1 ] 



(mo) 



G%\L,L)<Jr)t 



I \/Ti" } ) + /3 



(B21) 



for any positive numbers a and /3. Here, we used the relation S^oyT^ = S^a for any positive number a. By Eq. ( |B20t we 
also rewrite the quantity dB 1 31 > as 



°(",/3) - 



{[(J-Aa) 



-ij(M>)j* [A^I-A^)- 3 



(vO) 



[A^I-A^)- 1 ] 
G { ^(L,L) 



[(/-A^^wfVTM.fVTk* 5 



(B22) 



for any positive numbers a and /?. By the definition of \j Y« the real part of V T« is non-negative, so that we have 



ef^l y Ya^J + (3 j =0, as well as d(a + /3) = and d(a — /3) = 5(a — /3), for any positive numbers a and /3. By 
these facts and Eqs. (IB lib . (IB21b and (IB22I) we obtain 



r/.r 



S^S^ia - (3) 



-sm(aL)sm(/3L)^ ( [(I - A Q ) j ji 



-11 ("0) 



^(f-Ag)- 1 ] 

G { <?\l,L) 



ttD[ 

(no) 



(/0 



W( \/T%\ f Vy^ 



Cm) 



(B23) 



for any positive numbers a and /3. 



For the quantity D 



-ip] 



To ) we have 
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J- a 



T 



for T^'T^ > 
for T^T^ < 



( T[") - Tg ) 



2muj 



6 icS-p + —-{j i - V ) , 



(B24) 



for any numbers a and /3 except in the zero-measure region 



satisfying the condition I 'J = ±yT^ ^ 0, i.e. 

for almost all numbers a and (3. By introducing the quantity 



where we used Eqs. Q, dB2l > and dB19l ), the relation 

D(~ix,iy) = D(x,y), d(x) = 5(x), d(x 2 — y 2 ) = [S(x — 
\y\) + 5(x + \y\))/(2\y\) and S(— x) = 8{x) for any real num- 
ber x and y 0). By Eqs. ©, <|B3} and <|B8]) we further 
note 



-w\ \IT { ;\{^T^ ] 



2muj 



x [c('-P' + —(ji-v) 



(B25) 



_ „i(ti 2 -^ 2 +2m™/S)B 2 t/(2m) 



ir5[a 2 - (3 2 + —^- n) - ixi a 2 - f + -—i 



and by Eqs. $U%, dm! . (IB231 (lB24b . dBBt and [A^] * = A (_n) , we obtain 



+ 00 



dx [$ a {x,t)]*<f>p(x,t) 



= %-/3)-^sin(cd,)sin(/3Z) £ { [(/ - A*)" 1 ] ( " 0) }* 

n— — oo 
+ 00 +oc ( 

>< e e \c%-»\t)\^-^[{i-Hr i r ) 



— oo v=— oo 



. A (-(»-M))^-^(t) [(/-A, 3 ) 



-n(^o) 



(B26) 



(B27) 



for almost all positive numbers a and /3. 



X(a) as 



+00 



x{ P ) = y, x 



ipn 



(B28) 



X n = — I dpX(p)e- l P 



(B29) 



Now, we introduce the discretized Fourier transformation for any function X n of integer n, noting the formula 
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J_ dp exp(inp) = 27r<5„o for any integer n. By using 

this transformation and the formula J2n=-oo e%pn = 27r<5( j o) 
for any p E [— it, tt], the discretized Fourier transformation 

(x7y)(j>) of the convolution (X*Y)„ = Yl=-oo ^„- v Y u 
for any functions X n and Y n of n is represented as 



(X*Y)(p)=X(p)Y(p) 



(B30) 



for — 7r < p < 7r, as known as the convolution formula of the 
Fourier transformation. By using this property we have 



-oo +00 



E E ^ ) (t)A(-^»[(/-A,)- 1 ] 



("0) 



fl — — O0 U— — OO 



+00 +00 

E E Af-^'C^-^Wta-A,)- 1 ] 

/i— — 00 z;— — 00 

I 



(B31) 



because the discretized Fourier transformation of the left-hand 
side of Eq. ( 1B3U is equal to the one of the right-hand side of 
Eq. dB3TT >. By inserting ( tB3Tb into Eq. (lB27l > we obtain Eq. 
( TOl i for almost all positive numbers a and /?. 



Appendix C: Survival Probability under the initial conditions 
C0DandiI& 



In this appendix we give an explicit form of the survival 
probability P(t) under the initial conditions (fl~8T > and (120) . 
We note the integrals 

/•£ ^ 

/ da; sin(ax) sin(/3a;) = — S(aL, /3L), (CI) 
Jo 2 



cfo G^(o;,L)sin(^a;) 

. / T (nl , 



2VTi n) 



51 y/r£>L,pL 



dx 



G<£\x,L)\ CP(x,L) = 



Le 



T 



5( ( VT^) L^/TpL 



where <S(X, F) is the function of X and Y, defined by 



S(X,Y) = I 
2z 



i(-X--y) _ -i(X-Y) p i(X+Y) _ -i(X+Y) 



X -Y 



X + Y 



(C2) 



(C3) 



(C4) 



and satisfies the relation S(X, Y) = S(Y, X). By using Eqs. ©, ( fT2b . dCU and dC2| ), for the initial state given by Eqs. ( fl~8T > and 
(1201 the coefficient (TOT l is represented as 



Ai 



<S(fcL, cttt) + ^ K{k,n;0)S ( Tjj, n) , cttt 



(C5) 



where we define the function 7£(fc, n; t) of fc, n and t by 

[A fc (J - Afe)" 1 ] ( ™ 0) sin(fcL) .(^i-n^t) 



(C6) 



By Eqs. ©, (Hit . (fl4T >. dCTT l, dC2l , (IC3l . ( tC5l and dC6>, the survival probability (O is represented as 

r+oo r+co 

P(t)= da dp (A a )* A p B aP (t) (C7) 

Jo Jo 
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where B a p(t) is defined by 



+00 



B n 



-\s(aL,(3L) + 



fl— — OO 



* + OG 



y— — OQ 



fl— — OO V — — OO 



(C8) 



Therefore, the calculation of the survival probability P(t) B a p(t), and the double integrals ( IC7b - 
is attributed into the calculations of the quantities Ak and 



[1] G. Gamow, Z. Physik 51, 204 (1928). 

[2] R. W. Gumey and E. U. Condon, Phys. Rev. 33, 127 (1929). 

[3] H. J. Kupka, Transitions in Molecular Systems (Wiley-VCH, 
Weinheim, 2010). 

[4] D. J. Tannor, Introduction to quantum mechanics; A time- 
dependent perspective, (University Science Books, 2007). 

[5] R. G. Winter, Phys. Rev. 123, 1503 (1961). 

[6] W. van Dijk and Y. Nogami, Phys. Rev. Lett. 83, 2867 (1999); 
Phys. Rev. C 65, 024608 (2002). 

[7] G. Garcfa-Calderon, I. Maldonado, and J. Villavicencio, Phys. 
Rev. A 76, 012103 (2007). 

[8] I. V. Zozoulenko and T. Blomquist, Phys. Rev. B 67, 085320 
(2003). 

[9] G. Casati, G. Maspero, and D. L. Shepelyansky, Phys. Rev. E 

56, R6233 (1997). 
[10] K. M. Frahm, Phys. Rev. E 56, R6237 (1997). 
[11] D. V. Savin and V. V. Sokolov, Phys. Rev. E 56, R4911 (1997); 

Phys. Rev. E 57, 7367 (1998) [errata]. 
[12] F. Barra and P. Gaspard, Phys. Rev. E 65, 016205 (2001). 
[13] M. Puhlmann, H. Schanz, T. Kottos and T. Geisel, Europhys. 

Lett. 69, 313 (2005). 
[14] T. Taniguchi and S. Sawada, Phys. Rev. A 84, 062707 (201 1). 
[15] A. del Campo, F. Delgado, G. Garcfa-Calderon, J. G. Muga, 

and M. G. Raizen, Phys. Rev. A 74, 013605 (2006). 
[16] T. Taniguchi and S. Sawada, Phys. Rev. E 83, 026208 (2011). 
[17] A. del Campo, Phys. Rev. A 84, 012113 (2011). 
[18] D. Waltner, Semiclassical approach to mesoscopic sys- 



tems: Classical trajectory correlations and wave interference, 

(Springer- Verlag, Berlin, 2012). 
[19] L. A. Bunimovich and C. P. Dettmann, Europhys. Lett. 80, 

40001 (2007). 
[20] H. Massmann, Am. J. Phys. 53, 679 (1985). 
[21] D. Onley and A. Kumar, Am. J. Phys. 59, 562 (1991). 
[22] F. -M. Dittes, H. L. Harney, and A. Muller, Phys. Rev. A 45, 

701 (1992); H. L. Harney, F. -M. Dittes, and A. Muller, Ann. 

Phys. (N.Y.) 220, 159 (1992). 
[23] M. Miyamoto, Phys. Rev. A 68, 022702 (2003). 
[24] J. H. Shirley, Phys. Rev. 138, B979 (1965). 
[25] L. E. Reichl, The transition to chaos: Conser\>ative classical 

systems and quantum manifestations, (Springer- Verlag, New 

York, 2004). 

[26] M. Razavy, Quantum theory of tunneling, (World Scientific, 
Singapore, 2003). 

[27] Y. V. Nazarov and Y. M. Blanter, Quantum transport: Introduc- 
tion to nanoscience, (Cambridge University Press, Cambridge, 
2009). 

[28] M. V. Moskalets, Scattering matrix approach to non-stationary 
quantum transport, (Imperial College Press, London, 2012). 

[29] Our numerical calculations confirm the condition P(0) = 1 
for the initial conditions d 1 8b and l !20t . This is a supportive ev- 
idence that the initial state described by Eqs. j I St and J20b is 
expanded by the scattering states {&k(x, t)}k, namely the con- 
dition dT4t . 



